Quantum spatial correlations in high-gain parametric down-conversion measured by 

means of a CCD camera 
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We consider travelling-wave parametric down-conversion in the high-gain regime and present the 
experimental demonstration of the quantum character of the spatial fluctuations in the system. 
In addition to showing the presence of sub-shot noise fluctuations in the intensity difference, we 
demonstrate that the peak value of the normalized spatial correlations between signal and idler lies 
well above the line marking the boundary between the classical and the quantum domain. This effect 
is equivalent to the apparent violation of the Cauchy-Schwartz inequality, predicted by some of us 
years ago, which represents a spatial analogue of photon antibunching in time. Finally, we analyse 
numerically the transition from the quantum to the classical regime when the gain is increased 
and we emphasize the role of the inaccuracy in the determination of the symmetry center of the 
signal/idler pattern in the far-field plane. 



I. INTRODUCTION 



Quantum noise reduction has been the subject of ex- 
tensive research in the past. Although most of the 
experimental effort has been concentrated on the gen- 
eration and applications of quadrature-squeezed light, 
amplitude-squeezed light and intensity correlated twin 
beams have also been produced. The possibility of quan- 
tum correlations between light beams has played a ma- 
jor role in the recent development of quantum optics. In 
particular the spatial aspects of correlations of quantum 
optical fluctuations have been the object of several stud- 
ies (see for example Q, 0, S El IE @) an d have several 
new and promising applications such as quantum holog- 
raphy pi, the quantum teleportation of optical images 
P| , and the measurements of small displacements beyond 
the Rayleigh limit . An overview of this relatively new 
branch of quantum optics, defined as quantum imaging, 
can be found in ^Jj . 

The process of parametric down-conversion (PDC) is 
particularly suitable for the study of spatial effects, be- 
cause of its large emission bandwidth in the spatial fre- 
quency domain |Tlj . There is now a large literature 
on spatial effects in the low-gain regime, where photon 
pairs are detected via coincidence counting 0]. Never- 
theless, to date spatial correlation measurements in this 
regime have not evidenced any relevant quantum effect 
[l2l f]~3j - On the other side, measurements performed 
in the medium-gain regime (pump power < 1 MW) ev- 
idenced the twin beam character of the PDC emission 
[lij . i.e. a sub-shot-noise correlation between the pho- 
ton numbers of the whole signal and idler beams. Re- 
cent theoretical investigations done for an arbitrary gain 
[Tol Il6[ have predicted multi-mode spatial correlations 
below shot-noise between several portions of the signal 
and idler emission cones that correspond to phase conju- 
gate modes. There is a minimum size of the modes that 
have to be detected in order to observe a quantum corre- 
lation, which we shall refer to as the coherence area. This 
is determined by the conditional uncertainty in the direc- 



tions of propagation of the twin photons, and is roughly 
given by the inverse of the near-field gain area |15| . 

In this paper we present a description of the far 
field detection of the PDC radiation emitted by a j3- 
barium borate (BBO) non-linear crystal pumped by a 
low-repetition rate (2 Hz) pulsed high-power laser (1GW- 
lps). The detection is performed by means of a high 
quantum efficie ncy (77 rj 89%) scientific CCD camera. In 
a recent letter [13 we presented measurements showing 
the existence of quantum spatial correlation between the 
signal and idler beams. The aim of this paper is, on the 
one side, to describe in details several features of the ex- 
periment (such as the asymmetry of the PDC beam dis- 
tribution, the CCD diagnostic) which were not presented 
in 17]. On the other side, we shall provide a careful nu- 
merical analysis and interpretation of the experimental 
data, which (i) will enlighten a different aspect of the 
quantum noise correlation -namely an apparent violation 
of the Cauchy-Schwarz inequality-, and (ii) will discuss 
the role of the finite size of the detection pixel and of the 
inaccuracy in the determination of the symmetry center 
in the transition from the quantum to the classical regime 
of correlation. 

The use of the pulsed high-power laser enables us to 
tune the PDC to the high-gain regime while keeping a 
large pump beam size (of the order of ~1 mm). Thanks 
to the huge number of radiation transverse modes, we 
can concentrate on a portion of the parametric fluores- 
cence close to the collinear direction and within a narrow 
frequency bandwidth around degeneracy. This portion 
still contains a large (>1000) number of pairs of sig- 
nal/idler correlated phase-conjugate modes, propagating 
at symmetrical directions with respect to the pump in 
order to fulfill the phase-matching constraints. In the 
far field, where the measurement is performed, the cou- 
ples of modes correspond to pairs of symmetrical spots, 
which can be considered as independent and equivalent 
spatial replica of the same quantum system. Thanks to 
the very large number of these, the statistical ensem- 
ble averaging necessary for the quantum measurement 
can be solely done over the spatial replicas for each, sin- 
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gle, pump-laser pulse. Thus, differently from the exper- 
iment in [T^j . where the statistics was performed over 
different temporal replica of the system, here no tempo- 
ral averages over successive laser shots are considered. 
In our experiment the single-shot measurements reveal 
sub-shot-noise spatial correlations for a PDC gain corre- 
sponding to the detection of up to ~ 100 photoelectrons 
per mode |l7j . In connection with the work of |l8j , where 
the correlation between signal/idler phase-conjugate spa- 
tial modes was measured using a band-pass optical para- 
metric amplifier, a transition from the quantum to the 
classical regime with increasing gain was observed. Here 
we observe a similar transition, which is attributed to 
two main effects, one being the narrowing of the near- 
field gain profile that occurs at very high gain in presence 
of a bell-shaped pump beam. In fact this gain narrow- 
ing effect implies a broadening of the far-field coherence 
area, whose transverse size is measured experimentally 
from the width of the spatial intensity correlation func- 
tion evaluated between the symmetrical signal and idler 
regions under consideration. The other effect is the un- 
avoidable inaccuracy in the experimental determination 
of the symmetry center of the signal/idler pattern in the 
detection plane, and the role of the CCD pixel size in 
relation to the coherence area dimension of the radia- 
tion. Our experimental results are in accordance with 
the predictions of a 3D quantum model used to describe 
the experimental system |15). 

The paper is organized as follows. In section 2 we 
present the characteristics of the far field PDC radiation 
in the pumping conditions of the experiment, by illustrat- 
ing the features of the signal and idler beams detected by 
the CCD camera at degeneracy. In section 3 we carefully 
describe the experimental set-up used for the spatial cor- 
relation measurements and illustrate the typical single- 
shot images recorded by the CCD for the evaluation of 
the correlation, and in section 4 we present the experi- 
mental results. Our attention is focused on the analysis of 
the spatial correlation between signal and idler, to show 
that the height of the correlation peak demonstrates the 
quantum nature of spatial correlation and implies the 
apparent violation of a Cauchy-Schwartz inequality, i.e. 
a spatial analogue of photon antibunching in time, as 
predicted in |l6t IT^. I20I l2*lj . The numerical simulations 
showing the expected transition from the quantum to the 
classical regime for increasing gain are presented in sec- 
tion 5, where we also highlight the important role of the 
resolution cell size of the CCD camera (pixel size) and 
thus, as mentioned above, of the inaccuracy in the ex- 
perimental determination of the center of symmetry of 
the signal and idler image recorded in the far-field plane. 
The conclusions are given in section 6. 



(b) : Pump pulses 

(352nm, lps) 

Teflon pinhole 

(200 mm 

diameter) 



Type II 

rp BBO 



(a) 



Type II 
Pump BBO 
Pulses 



IK 



Mirror . IP 

WS3S2 ^ 5cm 



CCD 
Camera 



Rectangular 
aperture 

PBS 



/ A 



^4 



f=10cm 



CCD 
Camera 



FIG. 1: (a) Scheme of the diagnostics for the far-field de- 
tection of the degenerate signal and idler ring-type modes, 
(b) Detailed scheme of the experimental set-up used for the 
spatial correlation measurements. The third harmonic of the 
Nd:Glass laser is used to pump at 352 nm the BBO crystal 
which is cut for degeneracy at 704nm (9 = 49.05°, (f> — 0). 



II. DETECTION OF THE SPATIAL FEATURES 
OF THE FAR-FIELD PDC RADIATION BY 
MEANS OF THE CCD 



Before the quantitative investigation of the existence 
of spatial correlations between signal and idler beams, we 
perform a preliminary characterization of the generated 
parametric radiation. The type II 5x7x4mm 3 BBO non- 
linear crystal, operated in the regime of parametric am- 
plification of the vacuum-state fluctuations, is pumped 
by the third harmonic (352 nm) of a lps pulse from a 
chirped-pulsc amplified Nd:glass laser (TWINKLE, Light 
Conversion Ltd.). The input and output facets of the 
crystal are anti-reflection coated at 352 nm and 704 nm, 
respectively. The pump beam (vertically polarized (e)) 
is spatially filtered and collimated to a beam waist char- 
acterized by a full width at half maximum (FWHM) of 
approximately 1 mm at the crystal input facet, although 
the fluorescence characterization has also been performed 
for smaller beam sizes. The energy of the 352 nm pump 
pulse can be continuously tuned in the range 0.1-0.4 mJ 
by means of suitable attenuating filters and by changing 
the energy of the 1055nm pump laser pulse, allowing to 
have a gain G (representing the intensity amplification 
factor) in the range 10 < G < 10 3 . The parametric fluo- 
rescence of the horizontally polarized signal (o) and verti- 
cally polarized (e) idler modes is emitted over two cones, 
whose apertures depend on the specific wavelengths (see, 
e.g., lH HI). The BBO crystal (6 = 49.05°, = 0) is 
oriented in order to generate signal and idler radiation 
cones tangent to the collinear direction at the degener- 
ate wavelength uj s — uji — uj p /2 (s, i and p referring to 
signal, idler and pump respectively). 
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FIG. 2: Experimental far-field images of the degenerate sig- 
nal (left ring) and idler (right ring) beams recorded in single 
shot by the CCD placed in the focal place of a single lens 
(f=50 cm), with pump intensity (a) I » 30 GW/cm 2 , and (b) 
I « 50 GW/cm 2 , and FWHM pump beam size of 1 mm (a) 
and 0.4 mm (b) respectively. 



A simple far-field detection set-up is initially mounted 
as shown in Fig. Qfa). A deep-depletion back illumi- 
nated charged coupled device (CCD) camera [24| (Roper 
Scientific, NTE/CCD-400EHRBG1, with quantum effi- 
ciency ?yw 89% at 704 nm) triggered by a pulse from the 
laser system, is placed in the focal plane of a single large- 
diameter lens (f=5 cm), which collects at a distance f the 
far-field PDC radiation emitted by the BBO. The CCD 
detection array has 1340 x 400 pixels, with a pixel size of 
20/im x 20/xm. The pump-frequency contribution is re- 
moved by using a normal incidence high-reflectivity (HR) 
mirror coated for 352 nm placed after the BBO. By using 
a 10-nm broad interferential filter (IF), centered at 704 
nm, we are able to visualize the degenerate signal and 
idler far-field beams emitted in the parametric process. 
It is worth pointing out that without any spectral filter- 
ing, emission occurs on a very wide range of wavelengths 
and emission angles (see, e.g. j2^,|23j). Typical far-field 
images recorded at degeneracy in a single shot (for lps 
pump pulse) are shown in Fig.|2fa) and Fig.^b) for two 
different values of the pump intensity and, in the partic- 
ular cases illustrated, for two different pump beam sizes. 
The ring-shaped angular distribution is determined by 
the phase-matching conditions j^], and the rings at de- 
generacy are characterized by an angular width of about 
8° each. Note that with this set-up the two rings, which 
are emitted along the vertical direction, are recorded by 
rotating the CCD by 90° in order to fit the entire ring 
pattern inside the rectangular chip. 

We draw now the attention on the speckle-pattern as- 
pects of these images, observing that the dimensions of 



the speckles or spots (corresponding to the transverse 
modes in the far field) considerably change as a func- 
tion of the pump beam waist and in general of the pump 
beam intensity, and therefore as a function of the gain. 
The enlargement of the coherence area of the PDC radi- 
ation in the far field for increasing gain is evident from 
Fig. [21 and is confirmed by the theory and the numerical 
calculations. This effect, which will be extensively ana- 
lyzed in the following in sections 5 and 6, turns out to be 
very important for the interpretation of the experimental 
results of the correlation measurements. 

A closer inspection of Fig. [21 reveals that the photon 
distribution on the rings is asymmetrical: there are more 
photons along the outer edges of the rings than in the 
collinear region where the rings overlap. Moreover, this 
asymmetry seems to become stronger in Fig. Htb) which 
has a narrower pump waist than Fig. |2Ja) . An explana- 
tion of this phenomenon could be found by considering 
that for a small pump waist hot spots may appear in the 
emission spectrum of PDC. This is due to an additional 
phase-matching condition appearing when the Gaussian 
shape of the pump is taken into account |26j ; while always 
present when the pump field is Gaussian, it only becomes 
relevant when the pump waist is quite small. Instead, for 
larger pump waists the usual phase-matching conditions 
dominate. The hot spots are predicted to occur along the 
walk-off direction (in our case in the x direction in the ref- 
erence frame of Fig. [2J while perpendicular to the walk- 
off direction less photons should appear. To investigate 
if the hot spot phenomenon could explain the observed 
asymmetry we present in Fig.|3^a)-(b) numerical simula- 
tions for the same parameters as in Fig. [21 (for details on 
the numerics, see Ref. [l5j and section 5). Despite the 
hot spot centers being predicted to be located close to 
the region where phase matching occurs, no asymmetry 
in the photon distribution is observed in Fig. EHa)-(b). 
In order to see hot spots we had to reduce the pump 
waist further, as shown in Fig.^c): now hot spots along 
the walk off direction are evident. Figure Eld) shows the 
signal and idler fields plotted separately after averaging 
over 20 pump pulses, and a weak asymmetry between the 
central (collinear) region and the edges of the circles is 
seen. 

The numerics suggest that the hot spot phenomenon 
does occur and that it might lead to an asymmetry in the 
photon distribution. However, the pump waist needed to 
observe this is almost an order of magnitude lower than 
that used in the experiment. Therefore we checked if im- 
perfections in the pump beam reduce the effective waist, 
thus provoking the hot spot effect: simulations includ- 
ing a Gauss-Hermite pump shape with an overall large 
pump waist but with smaller spikes in the profile showed 
that hot spots could occur for larger pump waists than 
those found from a pure Gaussian beam. Thus a clean 
pump seems vital as to avoid hot spots from occurring, 
in particular considering that in the crystal Kerr nonlin- 
ear effects would amplify any spikes in the pump profile 
(an effect we have not accounted for in the numerics). 
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FIG. 3: Numerical far-field images calculated in a single shot, 
(a) and (b) are characterized by the same pump intensities as 
Fig. [5f a) and (b) respectively, and the pump intensity FWHM 
is (a) 1 mm and (b) 0.4 mm. (c) has pump intensity / ~ 
100 GW/cm 2 and a pump intensity FWHM of 0.12 mm. (d) 
shows an average over 20 shots of the signal and idler fields 
plotted separately for the same numerical conditions of (c). 
A 16 nm IF was used centered on 704 nm. From the theory 
in Ref. (26^1 . the predicted hot spot centers for signal (idler) 
are marked with squares (plus). 9 — 48.90° was used. 



In other words small-scale focusing of the pump (i.e. fil- 
amentation) might occur or at least be initiated in the 
crystal due to input-beam imperfections and Kerr nonlin- 
ear response. This might cause the spectrum to behave as 
if the actual pump-beam size is smaller than its FWHM 
diameter. Nonetheless, we do not believe that hot spots 
alone can explain the observed asymmetry. One alterna- 
tive explanation might be that absorption in the crystal 
is saturating the collinear direction of phase-matching as 
shown for a type I optical parametric amplifier by Dou 
et al. 27]. However, this seems an unlikely explanation 
because the BBO crystal has a very low absorption, and 
also Ref. [271 reports that the saturation occurs for much 
longer crystals than our. Preliminary numerical simula- 
tions that take into account such a loss also did not show 
any qualitative difference with respect to the data shown 
in Fig. |3 while it remains to be investigated if nonlin- 
ear absorption could lead to more dramatic effects. Wc 
are also suspecting that the interference filter is creating 
some unwanted attenuation of the collinear direction, or 
distortion effects, and we are currently investigating this 
in further detail. 



III. EXPERIMENTAL SET-UP FOR SPATIAL 
CORRELATION MEASUREMENTS 

The existence of spatial correlation already appears 
from the symmetrical properties of the signal and idler 
patterns recorded experimentally and shown in Fig. [21 
However, to investigate and reveal the quantum char- 
acter of the correlation we use the experimental set-up 
illustrated in Fig.^b) where, with respect to Fig. ^a), 
a different diagnostics configuration is adopted. We con- 
sider a pump beam size of 1mm and we now select the 
fluorescence around the collinear direction by means of a 



5mm x 8mm aperture placed 15 cm from the output facet 
of the BBO. The radiation is then transmitted through a 
polarizing beam splitter (PBS) that separates the signal 
and idler beams. The aperture prevents beam clipping 
by the PBS and thereby reduces substantially scattered 
radiation. The beams are finally sent onto two separate 
regions of the high efficiency CCD, which is placed in the 
common focal plane of the two lenses (f=10 cm) used to 
image the signal and idler far fields. In contrast to the 
case of photon-counting experiments (which rely on co- 
incidence measurements), and also differently from what 
was done for the preliminary characterization of the flu- 
orescence, in this set-up the correlation measurements 
are performed without using any narrow-band IFs, since 
these unavoidably introduce relevant transmission losses 
reducing the visibility of sub-shot-noise correlations, and 
could also introduce distortion or even attenuation effects 
that, as commented in section 2, are still unclear. Here 
the pump-frcqucncy contribution is removed by using 
normal incidence (M5) and at 45° (M4) high-reflectivity 
(HR) mirrors coated for 352 nm placed before and after 
the PBS, respectively, and a low-band pass colour filter 
(90% transmission around 704 nm) placed in front of the 
CCD. Note that a second PBS (not shown in the figure) 
is placed in the arm of the (e) idler beam to remove the 
residual contribution of ordinary (o) radiation reflected 
by the first PBS (3%), and a further HR@352 nm mir- 
ror (M4) is placed in the signal arm at a suitable angle 
in order to balance the unequal transmission of radia- 
tion in the two arms. All the optical components (except 
the colour filter) have anti-reflection coatings at 704 nm. 
The estimated quantum efficiency of each detection line, 
which accounts for both the transmission losses and the 
detector efficiency, is 7y tot ~75%. 

It is worth pointing out that prior to the experiment 
for the detection of quantum spatial correlations, we per- 
formed a test of the capabilities of the scientific CCD 
camera to perform spatially resolved measurements of 
photon shot-noise. The CCD used for the diagnostic has 
in fact been calibrated pixel by pixel to compensate for 
the gain inhomogeneity of the pixels on the CCD chip, al- 
lowing the retrieval of Poissonian statistics of the spatial 
fluctuations of an uniform enlightening in the full range 
of the camera dynamics .25] . The procedure works effi- 
ciently for thermal as well as for laser sources, provided 
that the wavelength and the coherence properties of the 
source are chosen in order to avoid the formation of equal 
thickness fringes in the chip (etaloning effect). Calibra- 
tion has also allowed the comparison at the shot-noise 
level of images recorded at different places on the chip. 
At that stage, before the realization of this quantum cor- 
relation experiment, retrieving the shot-noise in the CCD 
full dynamic range using classical sources certainly paved 
the way to spatially resolved photon noise measurements 
at the sub-shot-noise level, and turned out to be a neces- 
sary step to demonstrate quantum properties of images 
by means of CCD diagnostic. 

Fig. 0Ja) shows a typical far-field image recorded in 
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FIG. 4: (a) Single-shot far field image recorded by the CCD 
for a pump intensity I ~ 30 GW/cm 2 . The spatial areas for 
statistics are delimited by the black boxes selected within the 
degenerate signal and idler modes, spatially localized from 
the single shot image recorded with the lOnm-broad IF (b). 
(c) Zoom of two symmetrical areas of the signal and idler far 
fields. 

a single shot in the experimental configuration of Fig. 
H^b), where a fairly broadband radiation (i.e., the one 
transmitted by the rectangular aperture) is acquired in 
the signal (left) and idler (right) branches. The selection 
of the desired temporal and angular bandwidth around 
degeneracy is made by temporarily inserting in front of 
the CCD a 10-nm wide IF around 704 nm, allowing us 
to locate the collinear degeneracy point (see Fig. H{b)). 
The data analysis is limited within two rectangular boxes 
(black frames in Fig0fa)) corresponding to an angular 
bandwidth of 20 mrad x 8 mrad and to a temporal band- 
width smaller than 10 nm. The selected regions contain 
4000 pixels each. Since the aim of this work is to investi- 
gate pixel-pair correlation, and since the size of the CCD 
pixel approximately corresponds to the physical size of 
a replica (coherence area), the ensemble is large enough 
to perform the desired statistics. A zoom of the selected 
areas is presented in Fig^fc), where the rather spectac- 
ular symmetry of the intensity distribution in the signal 
and idler branches shows the twin-beam character of the 
phase-conjugate modes. 

IV. DETECTION OF QUANTUM SPATIAL 
CORRELATION: SPATIAL ANALOGUE OF 
PHOTON ANTIBUNCHING IN TIME 

Each of the signal and idler far-field pattern taken 
separately, looks like a speckle pattern produced by a 
pseudo-thermal source, such as for instance a ground 



glass illuminated by a laser beam. When this thermal 
light is splitted by a macroscopic device as a beam- 
splitter, the two resulting beams show a high level of 
spatial correlation, which is however limited by shot-noise 
|28ll29ll30| . The spatial correlation of the signal and idler 
beams generated by PDC is instead of microscopic ori- 
gin, and is not limited by shot-noise. The aim of this 
experiment is to show the sub-shot-noise nature of the 
spatial correlation of the PDC beams. We are first inter- 
ested in the symmetrical pixel-pair correlation, which is 
evaluated experimentally by measuring the variance u 2 _ i 
of the PDC photoelectrons (pe) difference n s — rii of the 
signal/idler pixel-pair versus the mean total number of 
down-converted pe of the pixel-pair. This variance is 

0"a-i = (K - n if) - ( n s - n,) 2 (1) 

where the averages are spatial averages performed over 
all the symmetrical pixel-pairs contained in the chosen 
regions. Each single shot of the laser provides a dif- 
ferent ensemble, characterized by its pixel-pair average 
pe number (n s +rii), in turn related to the paramet- 
ric gain. In the experiment, ensembles corresponding to 
different gains are obtained by varying the pump-pulse 
energy. We note that the read-out noise of the detec- 
tor, its dark current, and some unavoidable light scat- 
tered from the pump, signal and idler fields contribute 
with a non-negligible background noise to the process. 
This is taken into account by appl ying a standard cor- 
rection procedure (see for example |3l|), by subtracting 
the background fluctuations a 2 from the effectively mea- 
sured variance cr 2 s+b ^_, i+b ^ °f the total intensity differ- 
ence (signal+background)-(idler+background) obtaining 
a 2 _ i = <J 2 s+b \_/ i+b \ — 2a 2 . This background noise, having 
a standard deviation of 7 counts (±0.1 from shot to shot, 
estimated by repeating the measurement with the same 
pump-pulse energy) is measured in presence of pulse illu- 
mination over an area of the same size of the acquisition 
area and suitably displaced from the directly illuminated 
region. The validity of the data correction procedure 
and the shot-noise level (SNL) calibration are made by 
sending in the set-up (with no crystal) through the PBS, 
a coherent pulsed beam (@704nm) linearly polarized at 
45°: we have verified, for different laser energies, that 
the intensity difference fluctuations from the two coher- 
ent portions of beams recorded on the CCD lie at the 
SNL only if this background noise is taken into account 
and subtracted from the rough data acquired. 

FigEl shows the experimental results where each point 
is associated with a different laser shot. The data are 
normalized to the shot-noise level, and their statistical 
spread accounts for the background correction. Although 
the noise on the individual signal and idler beams is found 
to be very high and much greater than their SNL (=(n s ) 
and (rii) respectively), we observe an evident sub-shot- 
noise pixel pair correlation up to gains characterized by 
(n s + 7ij)Ri 15 — 20. Since in that regime the observed 
transverse size of the coherence areas (i.e. of the modes) 
is about 2-4 pixels, this approximately corresponds to 
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FIG. 5: Intensity difference variance o 2 a _ i normalized to the 
SNL (n s +rii). Each point (white circle) corresponds to a 
single shot measurement where the spatial ensemble statis- 
tics has been performed over a 100 x 40 pixels region. The 
triangles (each one obtained by averaging the experimental 
points corresponding to a certain gain) and their linear fit il- 
lustrate the trend of the data in the region between (n s + rii) 
=8 and 20. 



100 pe per mode. The maximum level of noise reduction 
observed experimentally agrees with the theoretical limit 
(dotted line in Fig. |SJ) determined by the total losses of 
the system (~ 1 — r\ tot 

We can have an idea of the transverse size of the 
mode by looking at the standard two-dimensional cross- 
correlation degree 
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between all the angularly symmetrical signal and idler 
pixels contained within the black boxes (see Fig. |3J| . This 
can be plotted for instance as a function of the horizon- 
tal and vertical shifts of the recorded image on the CCD, 
keeping fixed the position of the boxes. In general | 7 < 1 
with 7=1 for perfect correlation. The transverse section 
of the correlation function 7 plotted as a function of the 
horizontal shift x (in pixel units) and obtained from four 
single-shot images corresponding to different gains is pre- 
sented in FigHJa)-(d). We can notice how the FWHM of 
the curves increases for increasing gain, clearly reflecting 
the increment of the speckles size (and thus of the trans- 
verse mode size) already observed in Fig|21 For instance 
the graph obtained from a single-shot image character- 
ized by (n s + n,)« 8 (see Fig^a)) reveals a transverse 
mode size of about 2 pixels, while in the last case plotted 
in the figure where (n s + n^tx 40 the FWHM of the cor- 
relation function reaches a dimension of almost 7 pixels. 
As expected, virtually perfect correlation (in our case we 
have peak values of 7 up to ~ 0.99) is obtained for per- 
fect determination (i.e. within one pixel) of the center of 
symmetry between the signal and the idler regions. Fig[7| 
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FIG. 6: Correlation degree profiles plotted for four different 
gain values. 



illustrates the trend (as a function of the gain) of the co- 
herence area transverse size evaluated from the FWHM 
of the correlation degree profiles. The experimental re- 
sults (here associated with six images, the same as in 
FigElwith in addition two recorded for higher gains) are 
compared with numerical data calculated for two differ- 
ent pump size conditions. Although the qualitative trend 
between experiment and numerics is in good agreement, 
we observe a quantitative discrepancy between the ex- 
perimental and numerical data values. This discrepancy 
is attributed cither to an overestimation of the pump 
beam size in the laboratory (measured at a different time 
from the data acquisition period -e.g. one month before-, 
and possibly measured for different conditions of the laser 
source, which is very sensitive to laboratory temperature 
fluctuations and very sensitive to the laser regenerative 
amplifier misalignment); or to a suspected narrowing of 
the third-harmonic pump profile for very high gains that 
could occur as a consequence of the nonlinear genera- 
tion process, since the pump energy increase is performed 
by increasing in fact the fundamental laser pulse energy. 
This is a point which is currently under further investi- 
gation and new pump beam characterizations are going 
to be performed. 

It is interesting to note that the quantum nature of the 
correlation can also be estimated from the peak value of 
the correlation degree. As a matter of fact, since 



= {(n s - rii - (n„) - (rii)) 2 ) = 
o\ + of - 2((n s m) - (n s )(rii)) 



(3) 



the sub-shot noise condition for the intensity difference 
variance 



can be rephrased in the form 
7 > 1 



(4) 



(5) 



- Exp. results (estimated FWHM of pump ~ 1 mm) 

- Num. results (FWHM of pump - 1 mm) 

- Num. results (FWHM of pump ~ 600 (xm) 
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FIG. 7: Experimental and numerical results for the transverse 
mode size evaluated from the FWHM of the correlation degree 
profiles for different gains. 



FIG. 8: Experimental correlation degree (triangles) measured 
from five signal/idler far-field images for different values of the 
gain. 



if we use @ and assume that 

(n s ) ~ (n t ) = (n s ,i) and a 2 s ~ of = a 2 st , (6) 

as also confirmed experimentally within a good approx- 
imation. We can rewrite (J5J further by taking into ac- 
count that, as is well known, the signal and idler beams 
taken alone display a thermal-like statistics. Therefore, 
a 2 ra of = a 2 sl ns (n s ,i)(l + (n Sti )/M) 0, where M is 
the degeneracy factor representing the number of spatial 
and temporal modes intercepted by the pixel detectors. 
In the condition of the experiment the pump duration 
is slightly longer than the PDC coherence time while the 
pixel area is smaller than the coherence area, so that M is 
expected to be only slightly larger than unity |3lJ . Using 
this, © becomes 



FigJSJ illustrates the trend of the peak values of 7 for dif- 
ferent gains (black triangles), extracted from the same 
images used for Fig0 The dashed-dotted curve corre- 
sponds to the quantum standard limit Ju m obtained by 
interpolation of the function jn m = 1 — (n s i) /of j, calcu- 
lated for different gains using the values for the mean and 
variance obtained from the experimental far-field pat- 
terns considered. The full line is the theoretical limit 
obtained from J7J by using M = 2.4 as a fitting pa- 
rameter. We thus observe a spatial quantum correla- 
tion whenever the value of 7 lies in the region above the 
theoretical quantum limit. This limit becomes stronger 
(e.g more demanding) as the gain increases. The ex- 
perimental correlation values obtained are, as expected, 
compatible with the trend of the data plotted in Fig|5J 
and highlight a quantum correlation region up to values 
of (n s + n.i) of about 20 pe (also corresponding to at least 
100 pe per mode). For instance the first three triangles 



on the left correspond to three images characterized by 
an intensity difference variance that is clearly below the 
SNL in Fig[3 while the other three triangles correspond 
to images that are characterized by an intensity differ- 
ence variance above the SNL. Similarly to the case of the 
points plotted in FigEl the background noise correction 
leads to a statistical spread of the measured data also 
for the correlation degree. In relation to this comment 
we mention then that the experimental points (black tri- 
angles) represented in Fig|§j have been extracted from 
images characterized, for a given gain, by the highest 
correlation found. 

Finally we note that if we multiply (JSJ by of ., with 
the help of J5J and © we obtain 

(nsii) - (n s )(rii) > <j 2 s i - (n s ,i) (8) 

i.e., by defining 

(Sn s Sni) ee (n s rii) - (n s )(n,), (9) 

and 

(: 6n% :) = a 2 s>i - (n M ), (10) 

where the symbol : : indicates the normal ordering, we 
find the following inequality condition 

(5n s Sn l )>(:5ni l :) = l(:5n 2 s :)(:5n 2 :)} 1 / 2 , (11) 

equivalent to the sub-shot-noise level condition ©, and 
which states that the cross-correlation between signal 
and idler is larger than the (normally ordered) self- 
correlation. This corresponds to an apparent violation 
of the Cauchy- Schwartz inequality. This effect, which 
was predicted in 0, |2(], 0] for the case of the optical 
parametric oscillator and then generalized ^(| to the case 
of the travelling-wave optical parametric amplifier, rep- 
resents a spatial analogue of the phenomenon of photon 
antibunching in time. 
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N (pixels) 

FIG. 9: Numerical calculation of a^-i (normalized to SNL) 
between symmetrical portions of signal and idler plotted as a 
function of the detection area represented by N x N binned 
pixels. Different curves correspond to different values of the 
gain characterized by the mean number of down-converted pe 
per pixel pair (n s + m). 



V. TRANSITION FROM QUANTUM TO 
CLASSICAL REGIME: NUMERICAL RESULTS 

In order to interpret the observed transition from quan- 
tum to classical regime in Fig- EI we present in Fig.[!|]the 
results of the numerical calculations. The full quantum 
model accounts for the two transverse and the temporal 
degrees of freedom with propagation along the crystal, for 
the angular and chromatic material dispersion up to the 
second order, for the finite spatial and temporal widths 
of the Gaussian pump pulse, and for the experimental 
quantum efficiency r) to t- It has been important to in- 
vestigate how these features affect the spatial quantum 
correlation phenomena initially predicted by the plane 
wave pump theory, also to identify the conditions under 
which they can be observed in the experiment. The quan- 
tum description of the PDC is performed by treating the 
pump as a known classical field which propagates linearly 
inside the crystal, while the down-converted fields are 
quantized. The quantum averages needed in the calcu- 
lations are evaluated through a stochastic method based 
on Wigner representation, as described in |l5j . 

Fig- EH presents <r'^_ i , normalized to the SNL, vs. the 
size of the detection area for different gains. Each point 
is the result of a statistics performed over one single 
laser-shot. The case N=l corresponds to the experiment. 
The simulations (data not shown) outline that, in spite 
of the fixed pump-beam diameter, the signal and idler 
beam diameters at the crystal output strongly depend 
on the gain and decrease when the latter increases. This 
can be easily interpreted when considering that the sig- 
nal and idler beam size maps not the pump-beam pro- 
file but the actual parametric amplification gain profile 



G(r) ~ cosh 2 [crA(r)L] [33j (L being the crystal length, A 
the pump field amplitude and a a parameter proportional 
to the setting characteristics), as long as filtering due to 
the limited spatial bandwidth does not take place [34| . 
On narrowing the size of the PDC beams, the coherence 
areas in the far field (i.e. the modes) increase their size, 
as straightforward consequence of the convolution theo- 
rem in Fourier analysis 22]. This is precisely the effect 
observed experimentally and revealed either directly by 
looking at the far-field patterns recorded by the CCD, 
or by analyzing the correlation degree function profile as 
discussed in section 4. Since revealing quantum correla- 
tions requires detection areas larger (or comparable) to 
the mode size (as also discussed in [l5|), it is necessary 
when increasing the gain to have larger detectors in or- 
der to obtain sub-shot-noise variance as shown in Fig. 

Note that Fig. [5] evidences the transition from quan- 
tum to classical regime in case of single-pixel detection 
(N=l) for a gain that is higher than in the experiment. 
Indeed, in the experiment, excess noise is observed al- 
ready for (n s + rii) >20, which we attribute first to the 
effect of residual scattered light whose contribution grows 
linearly with the radiation fluence and is thus expected 
to overcome the shot noise at large pumping, and second 
to the uncertainty in the determination of the symmetry 
center of the signal and idler image portions. 

In fact although the increase of the coherence area with 
the gain leads to a deterioration of the signal-idler corre- 
lation, numerical calculations have shown that this phe- 
nomenon alone is not sufficient to explain the steepness 
of the slope given by the ratio a 2 _ i = d 1 s _ i l (n s + n s ) as 
a function of the shot-noise (n s + n s ) displayed in the 
plot of Fig. [S] As we shall now show, another impor- 
tant feature which contributes to this behaviour is re- 
lated to the inaccuracy in the determination of the cen- 
ter of symmetry of the far-field pattern in the detection 
plane. In our experiment the typical error Aa; symm be- 
tween the selected and the actual center of symmetry is 
on the order of half the size of the pixels of the CCD, 
that is about 10 /im. Ax symm can therefore be a sig- 
nificant fraction of the far-field coherence length, x co h 
(transverse size of the spatial mode), depending on the 
gain values considered in the experiment (here x co h is 
typically on the order of 50 /im). We are therefore in an 
intermediate situation between the ideal condition where 
A:c symm <C x co h ~ for which the theory predicts sub-shot 
noise correlation - and the opposite limit Ax symm 3> x co h 
in which the selected pixel pairs used to perform the 
statistics are completely uncorrelated. In this latter case 
the variance of n s — rii is simply the sum of the vari- 
ance of n s and rii (since (n s ni) = (n s )(rii)), that is 
crf_j = a 2 + of. As mentioned before, since the sig- 
nal and idler beams taken alone display thermal statis- 
tics, a 2 of = a 2 sl w (n s>i )(l + (n Syi )/M), having as- 
sumed (ri s ) (rii) = ( n s,i)- The saturation value of 
of_j in the Ax symm 3> x co h limit is therefore given by 
^sat = i°l + o-< )/(«- + ni) « 1 + (n s ,i)/M. 

The general behaviour for arbitrary values of the 
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FIG. 10: Plot of = a 2 a _ i j '{n s + in) as a function of 

Aa; symm for different value of the gain. For Aa; symm > x coh 
o'l-i saturates to a value on the order of 1 + (n s + rii}/ (2M), 
from which M ~ 3 is deduced. As a raw approximation the 
transition from sub-shot-noise to above shot-noise arises for 
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FIG. 11: Plot of al_i = v 2 a -i/{n s + m) as a function of 
{n 3 + rii) for different value of Ax symm . The slope of the 
curves is on the order of (2M)~ A% S ymm/%coh 



parameter is illustrated in Fig.^U] where a^_ i is 
obtained from numerical simulations by varying Aa; symm 
in the cc-axis (walk-off) direction for different gain val- 



ues. As Ax symm becomes larger than x co h the data 
saturate, and from these saturation values we estimate 
the degeneracy factor M rs 3 by comparing them to 
the predicted values of_ isat = 1 + {n s ^)/M. This M 
value is quite close to the one estimated experimentally 
(M = 2.4) in the previous section. From the behaviour 
of these curves it can be inferred that the transition from 
sub-shot-noise to above-shot-noise correlations occurs ap- 
proximately when Ax symm /x co h exceeds a critical value 

«l/^,sat«(l + K,,)/^)- 1 . 

It also follows from Fig. ^] that o^_ j increases almost 
linearly with {n s + rii) with a slope which is on the order 
of (2M)" 1 Ax symm /j; co ( 1 . The plot of Fig. HP confirms 
this result, showing how rapidly the sub-shot-noise cor- 



relation is lost as Axgymm becomes a significant fraction 
of x co h- We note that in the range of (n s +rii) considered 
in those simulations (which are those of the experiment) 
the effect of the broadening of the coherence area due to 
the increasing gain is negligible (see also the numerical 
resuls of Fig.[5]in the case N=l), as can be seen from the 
Ax symm = curve which is almost horizontal. For high 
values of the gain, the inaccuracy in the determination 
of the symmetry center becomes therefore particularly 
relevant and its contribution to the loss of correlation is 
expected to exceed the contribution due to broadening 
of the coherence area. Other contributions not included 
in the numerical model, such as light scattering from the 
environment, becomes also important as the number of 
PDC photons increases. 



VI. CONCLUSIONS 

In conclusion, we have presented experimental and 
numerical results demonstrating the quantum spatial 
features of the radiation generated in parametric-down- 
conversion in the high-gain regime and measured by 
means of a scientific CCD camera. We have first 
illustrated the speckle-like patterns of the fluorescence 
detected in particular at degeneracy, pointing out the 
effect of coherence area enlargement for increasing gain. 
The experimental investigation of the quantum aspects 
of the intensity correlations between signal and idler 
leads to the conclusion that twin beams of light gener- 
ated in PDC exhibit sub-shot-noise spatial correlation. 
This has been shown by measuring in the far-field an 
evident quantum noise reduction on the signal/idler 
intensity difference, and by having high peak correlation 
degree values lying above the standard quantum limit. 
The latter phenomenon corresponds to an apparent 
violation of the Cauchy-Schwartz inequality, i.e. to 
the spatial analogue of photon antibunching in time, 
which was predicted in |ll 13 E3, EH. A transition 
to above shot-noise correlation is observed as the gain 
increases. The theory and numerical simulations show 
how a quantum-to-classical transition is expected to 
occur because of a narrowing of the signal/idler beams 
with increased gain. This leads in turn to a larger 
far-field mode size and therefore also to the need of 
larger pixels to observe sub-shot-noise correlation [l5| . 
However numerical calculations have also shown that 
this phenomenon alone is not sufficient to explain 
the steepness of the slope of the transition observed 
experimentally. In addition to the effect of residual 
scattered light, whose contribution grows linearly with 
the radiation fluence and is thus expected to overcome 
the shot noise at large pumping, an important feature 
contributing to this behavior is found to be related to 
the inaccuracy in the determination of the center of 
symmetry of the signal/idler recorded pattern in the 
far-field plane. The importance of the results presented 
in this work lies in the fact that this is the first exper- 
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imental investigation of quantum spatial correlations 
in the high gain regime, where the huge number of 
transverse spatial modes is detected in single shot by 
means of a high-quantum-efficiency CCD. 
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